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Abstract 

Thermoelectric properties of a two-level molecule attached to the metal- 
lic electrodes are analyzed using the equation of motion technique within 
the Green's function formalism. Results show that the electrical con- 
ductance is strongly dependent on the electron and phonon temperatures 
and the electron-phonon coupling strength. In addition, it is observed 
that the thermal conductance peaks in the electron-hole symmetry points 
are vanished in the presence of the strong electron-phonon interaction. It 
is also found that the figure of merit is strongly suppressed in the strong 
electron-phonon interaction. The violation of the Wiedemann-Franz law 
is also observed coming from the Coulomb interactions. 

1 Introduction 

Transport phenomena through molecular devices have been attracted a lot of 
attention during recent two decades. Experimental and theoretical results show 
that the electron transport through the molecular transistors results in novel and 
interesting effects such as Coulomb and spin blockade effects [TJ [2j El H] , Kondo 
effect [SHI [7], negative differential conductance [H[9j[T0], and so on. In addition, 
it has been predicted that the molecular devices can be efficient for conversion 
of heat into electric energy [TTJ [T2]. Thermoelectric efficiency of a device is 
indicated by a dimensionless quantity called figure of merit, ZT = S 2 G e T/n, 
where G e is the electrical conductance and S is the thermopower. k is the 
thermal conductance which is the sum of the electronic contribution, K e , and 
the lattice thermal conductance, k p , and T denotes the operating temperature. 
In bulky samples and according to the Wiedemann-Franz Law, the ratio of the 
electrical conductance to the thermal conductance is related to the operating 
temperature. Research shows that the strong Coulomb repulsions, discreteness 
of energy levels, interference effects and etc. in nanostructures result in the 
violation of the Wiedemann-Franz law [T31 HH US] and as a consequence, the 
increase of ZT. For these reasons, the investigation of the thermoelectric prop- 



erties of the systems composed of a single or two quantum dots has attracted a 
lot of attention in recent years [Tll[Tg[ni[THl[Tll|2Dll2TJ|llll3IMIlS]. 

Electron-phonon interaction (EPI) is an interesting and important phenomenon 
in the molecular devices which can significantly affect the transport character- 
istics of the devices. The center of mass oscillation of the molecule [25], or 
thermally-induced acoustic phonons [27] can be the origin of the coupling be- 
tween the electronic degree of freedom and the vibrational degree of freedom. 
The influence of the EPI on the electronic characteristics of the molecular tran- 
sistors and the carbon nanotube quantum dots has been extensively studied 
using both rate equation approach [251 [2§1 1301 [33 131] and the Green's function 
formalism [331 [3H [3!3 13S1 [31] ■ However, the influence of the EPI on the thermo- 
electric properties of the molecular devices is an interesting issue needing more 
attention. Koch and co-workers |38j analyzed the thermopower in a single level 
molecule by means of the master equation. Furthermore, the effect of the EPI 
on the thermoelectric properties of a single-level quantum dot was studied in a 
few articles using the Green's function formalism [3J5] HUJ E3. • 

In this article, we consider a two-level molecule in which the electronic degree 
of freedom is coupled to the vibrational degree of freedom. Using the equation of 
motion technique within the Green's function formalism and using the polaronic 
transformation, the density of the states of the molecule (DOS) and the ther- 
moelectric properties of the device are analyzed. With respect to the fact that 
the phonon subsystem has a smaller heat capacity than the electron subsystem, 
it is completely probable an electron-phonon noncquilibrium thermodynamics 
is dominant so that the electron temperature, T e , and phonon temperature, T p , 
arc different. The difference can result in the appearance of the phonon absorp- 
tion sidebands in the DOS in the low electron temperatures which are not seen 
in the EP thermal equilibrium. The effects of the EP nonequilibrium and EP 
coupling strength on the thermoelectric properties of the system are analyzed 
in detail. In the next section, formalism is presented. Section 3 is devoted to 
the numerical results and in the end; some sentences are given as a summary. 

2 Model and formalism 

We consider a two-level molecule coupled to the normal metal electrodes. The 
Hamiltonian describing the system is given as 

aka i— l,2er ijcrcr' 

ua^a + Xn iAa> 1 + a] + £ K k *<ikA° + H - C -l 

ia otkrri 

where c akrT creates an electron with wave vector k, spin er, and energy £ a k<j in 
lead a. d| creates an electron in the ith level of the molecule with energy £j ff , 
while Ua stands for the on-site Coulomb repulsion and Uij (i ^ j) denotes the 
inter-level Coulomb repulsion. The energy levels of the molecule are tuned by a 
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gate voltage thus we set £.; CT = e° ia + V g so that Eq = e\ + A where A denotes the 
level spacing, a) (a) is the creation (annihilation) operator for the phonons with 
energy u> and A describes the electron-phonon coupling strength. The last term 
in the above equation describes the tunneling process between the electrodes 
and the molecule and V^ kcr is the coupling strength between the lead a and 
the ith energy level. V^ ktJ induces a tunneling rate from the molecule to the 
electrode a, T l aa = 27r^ feect p a \V^ k(J \ 2 where p a is the electronic density of the 
lead a. 

The EPI can be eliminated using polaronic transformation [55], H = e s He~ s , 
where S — exp(^,. (j \ni a [a' — a]). Thus, equation ([1} becomes 

H ^ ^ £akaC ak(J Caka ~t~ ^ ^ ^ia^ia ~t~ ^ ^ ^ UijTiifjTijtj (2) 
aka ia ijaa' 

akai 

where = ei a — A 2 /w, and Uij = — 2\ 2 /lj. The polaronic transformation 
results in the renormalization of the molecule levels and Coulomb repulsions. 
In addition, the tunneling amplitude is transformed into V^ ka = V^ ka X where 
X — exp(— ^[a^ — a]) is the phonon operator. In the following, we replace the 
operator X with its expectation value, < X >— exp(~^(N p + 1/2)) (55] where 
N p stands for the averaged number of phonons. This approximation is used in 
the localized polaron and is valid when the tunneling amplitude is smaller than 
the EPI, i.e. V a ka << A. Hence, the tunneling rate is renormalized according 
to T l aa = T l acr < X > 2 . < X > 2 just narrows the tunneling rate. 

In order to compute the electronic current, I e , and the heat flux, Iq, the 
Green function formalism is used so that the electronic current and heat current 
are expressed as [43l [44] 

_ P r°° pi pi 

I. = T E & ri £ lff to® - M^ ImG ^) ( 3a ) 

i poo pi 

T Q = IT £ / ^ ri L lr: ( £ ~ ~ /«(e)]/mGT„( e ) (3b) 

n ia J-oo 1 La + 1 i?.er 

where f a (s) = (1 + exp((e — /i Q ))//cT e ) _1 is the Fermi distribution function and 
/z Q , and T e are the chemical potential and electron temperature in the lead a. 
It is worth noting that equation ([3]) is computed with respect to the tunneling 
part of the Hamiltonian, equation (|TJ) , so the bare tunneling rates are used. 
However, the influence of the EPI on the electronic and heat currents is saved 
in the Green function, G\ a . -ImG r ia {e) is the DOS of the level i composed 
of the phononic and electronic parts which can be computed by evaluating the 
phonon part of the trace by Feynman disentangling technique |42j , and using the 
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Keldysh equations [33] for the lesser and greater Green functions, as follows [33J 



ImGUe)= £ ™ ""ft.- [rLUe + nu^maUe + nu) (4) 



L, 



■ni I pi 
a=L,.Rn=-oo Lo- T Act 

+ fL ff [l -/«(£- ™)]/mG^( £ - nw)], 
where L n , for nonzero temperatures, is given as 

L n = e -(V-) 2 (i+2^) ( ^±l ) ™/2 /n(2(A/w)2 J N {N )) (5 ) 

where I n (x) are the Bessel functions of the complex argument, and N p h = 
(exp(uj/kT p ) - I)- 1 is the averaged number of phonons and T p stands for the 
phonon temperature. Indeed, the thermal equilibrium between the electrons 
in the electrodes and phonons in the molecule can be broken due to smaller 
heat capacity of the phonons. G r icr (e) is the dressed retarded Green function 
governed by H. To compute G^ CT , we follow the procedure introduced by Chang 
and co-workers [45]. Therefore, G\ a is expressed as 

GL(e) = t^^Tr + E N r7 n } (6) 

f— * hj ia — llfc £ji a — Uu — Ilk 
k— u 

where Efo = e — E\ a + l/2(T l Lcr + r^ CT ), and, the summation is over all possible 
configurations in which the level j (j ^ i) is empty, po = 1+ < nj a rij- a > 
—(Nju+Nj-cr), singly, pi = N ja +Nj- a , or doubly, p 2 =< n ja nj- a >, occupied. 
Pk denotes the probability factor of a configuration expressed by one-particle 
Ni a =< dl a di a >g and two-particle < n, CT nj_ CT > i j occupation numbers. 11^ 
denotes the sum of Coulomb repulsions seen by an electron in level i due to 
other electrons in configuration k, in which the level j(j ^ i) is empty IIo = 0, 
singly, III = U12, or doubly, II 2 = 2C/i 2 occupied. 

In the linear response theory for the electronic current and heat flux (equa- 
tion (j3j), the electrical conductance is equal to G e = -y-Iu, thermopower is 

1 T 1 

S = -^ry 11 , and the thermal conductance is given by n e = ™[L22 — r 12 ]- 
The linear response coefficients are given by 46 

fa'lE/ ^ 'mGl,( c )[?M), (7b) 

" in J 1 La ~r 1 Ra u± e 



L ^ ??/ ^T^fft^^^ 2 ^^ (7c) 



For simulation purposes, we use half band width, D, as energy unit and set 
= D/50, and Uu = 1Un = D/10. We also take the renormalized tunneling 
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rate, f % a(7 as an input parameter and ignore the narrowing the width of the 
tunneling rate due to the EPI and assume the wave functions of the molecular 
levels are identical, i.e. T\ a — r^ CT . It has been shown that the EPI does not 
narrow the width of tunneling rate in some cases [UJ 135] . 
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Figure 1: (a) Electrical conductance as a function of electron temperature for 
A = (solid line), A = 1,T P = T e (dashed line), and T p = 5T e (dash-dotted 



line), (b)-(d) show DOS in the partsl-3, respectively. 
T l ary = 0.2w, and A = 2w is the level spacing. 



dm 

de 



is plotted in gray. 



3 Results and discussions 



Figure la describes the electrical conductance, G e , as a function of temperature 
in the absence and presence of the electron-phonon interaction. Results show 
that G e is uniformly increased by increase of temperature. However, behavior of 
it is strongly dependent on the electron and phonon temperatures and EPI. One 
can divide figure la to three parts. In part 1, low temperatures, kT e < 0.2oj, the 
electrical conductance resulting from the EP nonequilibrium is more. Figure lb 
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can help us understand the reason. In low temperatures the DOS of the molecule 
exhibits a vibrational absorption sideband (dash-dotted line) located near the 
chemical potential of the leads. Such peak is not observed in the elastic transport 
(solid line) or in the presence of the thermal equilibrium between electrons and 
phonons (dashed line). Note that there are no phonons in the molecule when 
phonons and electrons are in the thermal equilibrium because of kT p << u>. In 
region 2, mediated temperatures, 0.2 < kT e < uj, the electrical conductance of 
the elastic transport is more. However, EP nonequilibrium can enhance G e in 
comparison with the EP equilibrium. It is worth noting that G e resulted from 
the EPI is independent of the temperature difference between electrons and 
phonons for kT e = uj. In region 3, high temperatures, kT e > uj, G e approaches 
saturated values. The DOS of the molecule and (—/'(e)) are plotted in figures 
la-d for parts 1-3, respectively, /'(e) is a Lorentzian function whose center is 
located in \x and its width depends on the electron temperature. Therefore, 
/'(e) becomes wider with increase of the temperature so that more part of the 
DOS involves in the electron transport. 




Figure 2: Thermal conductance as a function of gate voltage. Inset shows the 
DOS when V g = -1/2(A + U x + 2U 12 ). T e = 0.3lu 
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The electrical conductance, K e , is plotted in figure 2 as a function of the 
gate voltage. In A = 0, the electrical conductance has seven peaks. Four 
peaks are located in resonance energies i.e. V g = e\, V g = —(A + U\ 2 ), 
v g = -(e° + Ui + U12), and V g = -(A + U 2 + 2U 12 ), in which the electrical con- 
ductance is maximum because electrons can easily tunnel to the molecule. Other 
peaks are located in electron-hole symmetry points, i.e. V g = — 1/2(A + U\ 2 ), 
V g = -1/2(A + Ui + 2U 12 ), and V g = -1/2(A + Ui + U 2 + 3t/i 2 ). In these 
points, thermopower is zero. In the electron-hole symmetry points, electrons 
and holes participate in the charge and energy transport with the same weight. 
Although they carry the charge in the opposite directions, they carry the en- 
ergy in the same direction resulting in the appearance of the peak in the thermal 
conductance spectrum. In addition, it is observed that the height of the peaks 
is strongly dependent on the electron population of the molecule, so that the 
peaks are lower in the voltages V g = e\, and V g = —(A + U 2 + 2Ui 2 ) which the 
molecule is, respectively, approximately empty or fully occupied. The peaks are 
higher when the molecule is partially occupied. The dependence of the DOS on 
the one- and two-electron populations gives rise to such behavior. In the pres- 
ence of the EPI, the height of the peaks is significantly reduced. Furthermore, 
peaks located in the electron-hole symmetry points are not well observed. Inset 
shows the DOS in V g = -1/2(A + Ui + 2U 12 ). It is clearly observed the DOS is 
completely symmetric around the chemical potential of the leads in the elastic 
transport, while the EPI disturbs the symmetry leading to the disappearing the 
peaks of the thermal conductance in the electron-hole symmetry points. How- 
ever, interaction of the EP nonequilibrium can slightly increase K e due to the 
existence of the phonon absorption side peaks. 

Figure 3 describes the behavior of the thermopower and figure of merit as a 
function of temperature in the absence and presence of the EPI. It is observed 
the thermopower is increasingly enhanced in the low temperatures, whereas 
it is reduced in the high temperatures. Because thermopower has the main 
contribution in the figure of merit, ZT follows the same behavior. Furthermore, 
results show that the thermopower is independent of the EPI if kT e = u. In 
low temperatures, holes participate in the generation of the current because 
£1 is below the resonance and electrons from the colder lead (right lead) can 
tunnel to the warmer one (left lead). Note that S in the unit of k/e is negative 
because of negative electron charge. With increase of temperature resulting in 
the excitation of the electrons toward above the chemical potential , electrons 
from the warmer lead can tunnel to the colder one using the level e 2 . It is so 
interesting that the sign of the thermopower changes faster in the presence of the 
EPI. It comes from the fact that the electrons and holes can more easily tunnel 
from one lead to the other by the phonon absorption processes. In addition, the 
magnitude of ZT is significantly reduced in the presence of the EPI. However, 
EP nonequilibrium can enhance ZT when 0.6w < kT e < lj because of the 
increase of the thermopower (figure 3a). ZT approaches the lower values in 
the high temperatures due to the inter-level and intra-level Coulomb repulsions 
(proximity effect). Such behavior was analyzed for a triple quantum dot system 
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Figure 3: (a) Thermopower and (b) figure of merit as a function of electron 
temperature. Parameters are the same as fig.l. 



by Kuo and co-workers |17j . 

The dependence of the figure of merit on the phonon temperature and EP 
coupling strength is analyzed in figure 4. In A = 0, ZT is independent of the 
phonon temperature so that its magnitude is constant and maximum. With 
increase of A the magnitude of ZT is significantly reduced by increase of T p . 
However, in weak As, A < OAlu, ZT is nearly constant but in the strong As, ZT 
is suppressed because of the reduction of the thermopower and increase of the 
thermal conductance. 

The Lorentz ration, L/Lq, where Lq = Tr 2 k 2 /3e 2 is plotted in figure 5 as a 
function of the molecule energy levels and in figure 6 as a function of temper- 
ature. The Coulomb interactions give rise to the violation of the Wiedemann- 
Franz law. It is observed that the Lorentz ratio oscillates as a function of the 
energy level. When the molecule is completely full or empty, the Lorentz ratio 
is one for elastic transport because the transport in this regime is dominant 
by higher order tunneling, while it is bigger than unity for the EP nonequilib- 
rium due to the phonon assisted tunneling. The Lorentz ratio is significantly 
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Figure 4: Color map of the figure of merit as a function of phonon temperature 
and EP coupling strength. kT e = 0.2w. 



increased in the Coulomb blockade regime. On the other hand, the dependence 
of the DOS on the one- and two- electron populations leads to the asymme- 
try of the Lorentz ratio. The EP nonequilibrium increases the Lorentz ratio 
near the chemical potential of the leads. The behavior of the Lorentz ratio as 
a function of temperature is different in the absence and presence of the EPI. 
In the presence of the EP nonequilibrium, the Lorentz ratio is lesser than one 
because the electrical conductance is increased due to the phonon absorption- 
assisted tunneling, while the Lorentz ratio is increased in the elastic transport. 
The reason of the difference was investigated in figure 1. In high temperatures, 
the Lorentz ratio decreases significantly because the electrical conductance is 
increased. Such behavior was previously reported for the single-level quantum 
dots [15]. 

4 Summary 

We analyze the thermoelectric properties of a two-level molecule by means of 
equation of motion technique within the Green function formalism. The influ- 
ence of the Coulomb interactions, electron-phonon coupling, and temperature 
on the figure of merit is examined. The electron temperature of the bulky elec- 
trodes can be different from the phonon temperature because of the smaller heat 
capacity of the phonons so that the electron-phonon nonequilibrium thermody- 



9 



Figure 5: Lorentz ratio as a function of gate voltage. Parameters are the same 
as fig. 2. 



namics may govern the system. The temperature difference results in the novel 
and interesting phenomena. It is observed that the electron-phonon interaction 
results in the reduction of the number of peaks in the thermal conductance 
due to the disturbance of the symmetry of the DOS near the chemical poten- 
tial of the leads. Furthermore, the Coulomb repulsions and electron-phonon 
interaction result in the violation of the Wiedemann-Franz law. 
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